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We study both analytically and numerically demagnetization via nucleation of the metastable 
phase in a two— dimensional nonequilibrium Ising ferromagnet at temperature T. Canonical equi- 
librium is dynamically impeded by a weak random perturbation which models homogeneous dis- 
order of undetermined source. We present a simple theoretical description, in perfect agreement 
with Monte Carlo simulations, assuming that the decay of the nonequilibrium metastable state is 
due, as in equilibrium, to the competition between the surface and the bulk. This suggests one to 
accept a nonequilibrium free-energy at a mesoscopic/cluster level, and it ensues a nonequilibrium 
surface tension with some peculiar low-T behavior. We illustrate the occurrence of intriguing 
nonequilibrium phenomena, including: («) stochastic resonance phenomena at low T which sta- 
bilize the metastable state as temperature increases; (m) reentrance of the limit of metastability 
under strong nonequilibrium conditions; and {in) noise— enhanced propagation of domain walls. 
The cooperative behavior of our system, which is associated to the interplay between thermal 
and nonequilibrium fluctuations, may also be understood in terms of a Langevin equation with 
additive and multiplicative noises. We also studied metastability in the case of open boundaries 
as it may correspond to a magnetic nanoparticle. We then observe the most irregular relaxation 
triggered by the additional surface randomness. In particular, at low T, the relaxation becomes 
discontinuous as occurring by way of scale-free avalanches, so that it resembles the type of relax- 
ation reported for many complex systems. We show that this results from the superposition of 
many demagnetization events, each with a well-defined scale which is determined by the curvature 
of the domain wall at which it originates. This is an example of (apparent) scale invariance in a 
nonequilibrium setting which is not to be associated with any familiar kind of criticality. 

KEY WORDS: metastable state, nonequilibrium, kinetic Ising model, nucleation, critical droplet, 
stochastic resonance, surface effects, disorder, scale-invariant noise, avalanches. 



I. INTRODUCTION 

Many different natural phenomena involve metastable states that, eventually, decay via nucleation. Some familiar 
examples appear in the flow of electrical current through resistors, relaxation in amorphous materials, glasses and 
gels, domain wall motion in hysteretic disordered magnets, granular media evolution, earthquake dynamics, protein 
conformations, or false vacuum states in quantum field theory. There is a great amount of information concerning 
these situations but a full microscopic theory of metastability and nucleation is elusive HI; 0; H; H 0; @; HSU)- To 
begin with, there are two main coupled difficulties. One is that metastability concerns dynamics (Isl;^). The systems 
of interest typically show a complex free-energy landscape with (many) local minima, which are metastable in the 
sense that they trap the system for a long time. One may imagine that, eventually, relaxation occurs when the system 
after long wandering finds a proper path between the minima. This results in a complicate coupling of dynamics and 
thermodynamics (jf Ot ). A second difficulty is that many systems of interest cannot reach thermal equilibrium even 
after relaxation. In general, they are open to the environment, which often induces currents of matter or energy, 
or they are subject to agents which impose opposing tendencies which typically break detailed balance (fill ). These 
■ perturbations result in a final steady state which cannot be described by a Gibbsian measure, i.e. a nonequilibrium 
. ^ , stationary state. Consequently, thermodynamics and ensemble statistical mechanics do not hold in these systems, 
' which is a serious drawback. 
H ' These difficulties make the field most suitable for simple-model analysis. Indeed, the (two-dimensional) kinetic Ising 



model has been the subject of many studies of metastability not only in the case of periodic boundaries ([Q; 12; 13; l4) 
but also for finite lattices with free boundaries (flU : [H; [TtI; ilB)- The latter try to capture some of the physics of 
demagnetization in very dense media where magnetic particle sizes typically range from mesoscopic down to atomic 
levels. If one keeps oneself away from specific models, metastable states are often treated in the literature as rare 
equilibrium states, at least for times much shorter that the relaxation time, and it has been shown that one may define 
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a metastable state in aproperly constrained (equilibrium) ensemble (@; H; @), and that most equilibrium concepts 
may easily be adapted S H H [1). 

In this paper we present a detailed study of metastability (and nucleation) in a nonequilibrium model. In order to 
deal with a simple microscopic model of metastability, we study a two-dimensional kinetic Ising system, as in previous 
studies. However, for the system to exhibit nonequilibrium behavior, time evolution is defined here as a superposition 
of the familiar thermal process at temperature T and a weak completely-random process. This competition is probably 
one of the simplest, both conceptually and operationally, ways of impeding equilibrium. Furthermore, one may argue 
that it captures some underlying disorder induced by random impurities or other causes which are unavoidable in 
actual samples. The specific origin for such dynamic randomness will vary with the situation considered. We mention 
that a similar mechanism has already been used to model the macroscopic consequences of rapidly-diffusing local 
defects (fill ) and quantum tunneling |l6l ) in magnetic materials, for instance. 

The weakest perturbation of that kind happens to modify essentially the canonical equilibrium. We observe, 
however, some structural similarity at the mesoscopic level which suggests that cluster dynamics is, as in equilibrium, 
the consequence of a competition between the surface and the bulk. Therefore, we postulate the existence of a 
nonequilibrium free-energy and an associated nonequilibrium surface tension. These functions are formally similar 
to the ones in equilibrium but behave peculiarly due to the different nature of the system states. In fact, two main 
predictions follow that are counter intuitive. One is that, at low T, the stability of the metastable state is enhanced by 
increasing temperature, as a consequence of resonance between thermal and nonequilibrium noises. The second main 
prediction is that the limit of metastability or pseudospinodal, which separates the metastable phase from the unstable 
one, exhibits reentrant behavior as a function of T for strong enough nonequilibrium conditions. Our predictions are 
nicely confirmed in related Monte Carlo (MC) simulations. We also show that the most important aspects of the 
model behavior are contained in a Langevin-like description which contains a multiplicative noise in addition to the 
more familiar additive one. 

The theory above concerns an infinite system, and the related computer simulations are, therefore, for lattices 
with periodic boundary conditions. However, we also studied other types of boundary conditions, including free 
borders. Previous work showed that free borders, which are most relevant when dealing with magnetic nanoparticles, 
for instance, determine importantly the particle demagnetization (fl7t [ish . We here confirm this and observe that, 
for free borders, but not for periodic boundary conditions, demagnetization occurs, via avalanches, through many 
different metastable-like configurations. These avalanches exhibit power-law size and lifetime distributions in a way 
that closely resembles the relaxation in many cases in nature. We show that this is a consequence in the model of 
the superposition of different events, each with a typical scale which is determined by the curvature of the interface 
at which the avalanche originates. 

The paper is organized as follows. We define the model in the next section. In section Hill we set the basis of our 
approach. In particular, ijllLAI concerns a mean-field approximation which allows for a discussion in i ^IILBI of the 
nonequilibrium surface tension. The next subsection introduces our ansatz for a nonequilibrium free-energy cost of 
excitations in the metastable phase, and we test the ensuing predictions against MC simulations. In i jlll.Dl we present 
a Langevin theory which is also a good description of the nonequilibrium situation of interest here. ijlH.El is devoted 
to study the limit of metastability. In section HVl we study the case of open boundaries. The statistics of the resulting 
avalanches are analyzed in ijlV.Al and i jlV.BI and i jlV.CI are devoted to discussions, respectively, of our observations 
and of its possible extension to interpret avalanche data from actual systems. Finally, section f|V] contains a brief 
conclusion. Some technical details are left for the appendices. 



II. THE MODEL 

Consider a two-dimensional square lattice of side L and, for the moment, periodic boundary conditions (or, even- 
tually, the thermodynamic limit condition L — )■ oo). We define a spin variable = ±1 at each node, i e [1, = i^]. 
Spins interact among them, and with an external magnetic field h, via the Ising Hamiltonian function 

N 

H{s)^- J2 -^.sj-hj^^^^ (1) 

\i-3\=l i=l 

where s = {si} and the first sum runs over all nearest-neighbors (NN) pairs. Time evolution proceeds by stochastic 
dynamics consisting of single-spin flips with the transition rate 

Lo{s^s')=p+{l-p)^iPAH^), (2) 

where < p < 1 is a parameter, s and s* stand, respectively, for the system configuration before and after flipping 
the spin at node i, ATCi is the energy cost of the flip, and p = 1/T. The (arbitrary) function in ^ is set ^'(a;) = 
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FIG. 1 (Color) Phase diagram for the nonequihbrium Ising ferromagnet of interest here. The sohd (red) hne corresponds to 
the critical temperature as obtained from the nonequihbrium surface tension in ijlll.BI The order disappears at T = for 
p > Pc = {y/2 — 1)^ ~ 0.17. Notice, however, the reentrant blob for pc < p < p^ ~ 0.186 and intermediate temperature. The 
dashed (blue) line is derived from a first-order (pair approximation) mean-field theory, namely, eq. ([8} in i jlll.AI 

e~^(l -I- e"^)"-'^, corresponding to the so-called Glauber rate, in our simulations. However, we shall also refer below 
to the case ^(a;) — min[l,e~^] which corresponds to the Metropolis rate. 

For any p, two different heat baths compete in One is at temperature T and operates with probability (1 — p), 
while the other induces completely-random {infinite temperature) spin flips with probability p. As a consequence of 
this competition, a noiiequilibrium steady state sets in asymptotically which cannot be characterized by a Gibbsian 
invariant measure (fill , [ig ). In particular, the rate ^ violates detailed balance for any p > 0. However, the system 
reduces for p = to the familiar Ising model with rate 4" which goes asymptotically to the equilibrium state for 
temperature T and Hamiltonian Ti.. This model is a particular case of a more general class of systems, characterized 
by a coupling to a number of different heat baths leading to nonequihbrium behavior, whose static properties have 
received considerable attention in recent years (fill). 

For zero magnetic field, the model exhibits an order-disorder continuous phase transition at critical temperature 
Tc{p) < Tc{p = 0) = Tons, where Tons — 2/ln(l + \/2) is the Onsager temperature. Fig. [T] shows the model 
phase diagram Tc{p) as obtained by two different approaches, as described in i jllLAI and ijlH.Bl respectively. In 
the former case, corresponding to a first-order (pair) mean-field approximation, the predicted equilibrium [p ~ 
0) critical temperature is the Bethe temperature Tbothc = 4/ln4 « 1.27Tons, and order disappears at T = for 

p > p i^^^^ = 5/32 « 0.156. A more accurate estimation follows from the nonequihbrium surface tension ctq in 
i jllLBI In this case, one recovers the correct equilibrium limit, and order is predicted to disappear at T = for 

p > pi""^ — {\/2 — 1)^ ~ 0.1716 (in agreement with MC results (fllf)). However, the anomalous behavior of (Tq that we 
describe in ijllLB I implies the emergence of an intermediate-T region for pc < p < p* ~ 0.18 625 where order sets in. 
This reentrant behavior of Tc{p), which is shown in Fig. [TJ is reminiscent of the one in the phase diagram for systems 
subject to multiplicative noise, as discussed below. 

III. NONEQUIHBRIUM METASTABILITY AND NUCLEATION 

In order to produce metastability in the above nonequihbrium setting, we initialize the system in a state with all 
spins up. Si = -1-1, i — 1 . . . N. For small h < and T < Tc{p), this configuration quickly relaxes to a metastable state 
with magnetization m — N~'^ Si^i > ^- There is, however, a tendency of the spins to line up in the direction of the 
field, which competes with the tendency to maintain the local order induced by the spins mutual interactions. The 
result for small \h\ is a metastable state of long lifetime, i.e., the system fiuctuates around the metastable minimum 
for a long time, though eventually decays toward the truly stable phase, which has m < 0. 

A. The Metastable Phase: Mean-Field Theory (Pair Approximation) 

The resulting metastable state is homogeneous at mesoscopic scales for times shorter than its relaxation time, as 
expected (0; 0), and we can exploit this feature. In particular, we shall assume that: (i) the system is completely 
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TABLE I Spin classes for the two-dimensional Ising model with periodic boundary conditions. The last column shows the 
energy cost of flipping the central spin at each class. 



homogeneous, and (ii) it can be described by the average density of up spins, p{+), and by NN (pair) correlations as 
captured by the density p{s, s') of NN spins. Under these assumptions, all the spins will behave in the same manner, 
and the evolution of a spin will depend exclusively on its current state and on its four NNs. Denoting x = p{+) and 
z = p{+, +), we may write p(— ) = 1 — x, p{+, — ) = p{+) — p{+, +) — x — z and p(— , — ) = 1 + z — 2x. In this 
approximation, the probability Q{s,n) of finding a spin in a state s — ±1 surrounded by n G [0,4] up NN spins is 



Q{s,n) = [ ^]p{s)p{+\srpHs)'- = 



p(+,5)>(-,s) 

n 



4~n 



Pisy 



(3) 



where p{±\s) = p{±,s)/p{s) are marginal densities. The pair of indexes (s,n) defines the spin class to which the 
given spin belongs to. For periodic boundary conditions, there are 10 different spin classes, as shown in Table L The 
cost A'H{s,n) of flipping any spin within a class is the same, so that the transition rate uj{s — > s*) = uj{s,n) in ([2|) 
only depends on s and n, the spin at a node i and the number of its up NN neighbors, respectively. 

This information may be used to write down rate equations for the two relevant observables in our mean field 
approximation, namely x and z. In particular, the change in x when a spin s with n up neighbors flips is Sx(s, n) = — s, 
whereas the change in z in the same process is 6z{s,n) = —sn/A. We thus obtain (jlli: i22): 



where 



G {x, z; n) 



dx 
di 

d£ 

dt 



Gi{x,z) EE ^G{x,z;n) , 



n=0 
4 



G2ix,z) = ^-G{x,z;n) , 



n=0 



(x - z)"(l + z - 2x)4-" , ^ z"(x-z)4-" , 



(4) 
(5) 

(6) 



Eqs. (IH)-® correspond to a mean-field Pair Approximation (jlll: 1231). This is a dynamic generalization of the first- 
order cluster variation method by Kikuchi ((2^ , and has been used to study other nonequilibrium systems with success 
([23I ). Our interest here is on the stationary solutions, G'i(a;st, Zgt) = G2(xst,Zst) = 0. Both stable and metastable 
solutions are locally stable against small perturbations, which requires the associated Lyapunov exponents to be 
negative. This yields the necessary and sufficient conditions 



(a,Gi)3,-l-(9.G2),, <0, 



(7) 



On the other hand, the condition (9a;Gi)st = signals for h ~ an incipient or marginal instability corresponding to 
an order-disorder phase transition with critical point (x^^, z^J = (1/2, 1/3). It follows that 



Tcip) = -4 



1 
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(8) 
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FIG. 2 (Color) Temperature dependence of the steady-state magnetization for h = and, from top to bottom, p = 0, 0.001, 
0.005, and 0.01. Lines are results from the mean-field (pair) approximation. Inset: Magnetization in the metastable phase vs. 
T for h = —0.1 and the same same values of p as in the main graph. 



This function is depicted in Fig. [T] (see the related discussion above in ^TT|). The stationary state (xstj-^st) may 
be obtained numerically from the nonlinear equations Gi{xst, Zst) = G2{xst, Zst) = subject to the local stability 
conditions in eq. ([7]). For h = 0, the up-down symmetry leads to pairs of locally-stable steady solutions, namely 
(xst, Zst) and (1 — Xgt, 1 + -^st — Sxgt)- The result is illustrated in Fig. [21 This also illustrates the expected agreement 
with Monte Carlo results at low and intermediate temperature for any p. In particular, this agreement is good for 
temperatures below 0.75 Tc{p). The fact that increasing p at fixed T decreases the magnetization reveals that the 
nonequilibrium perturbation tends to increase disorder, as expected. In this way we can identify the parameter 
p, which weights the competition between different heat baths leading to nonequilibrium behavior, as a source of 
nonequilibrium noise, as compared with thermal fluctuations induced by T. For small enough fields, the situation 
closely resembles the case h = 0, but the up-down symmetry is now broken, and locally-stable steady states with 
magnetization opposite to the applied field are metastable, see inset to Fig. [51 



B. Nonequilibrium Surface Tension 

The above mean-field approximation neglects fluctuations, so that it cannot account for the relaxation of the 
metastable phase. One expects that, for long enough times, fluctuations will allow for a decay toward the stable 
phase. Direct inspection of escape configurations as the ones in Fig. [3l indicates that the metastable-stable transition 
is a highly inhomogeneous process triggered by large stable-phase clusters, and a detailed comparison (not shown) 
suggests that the nonequilibrium and equilibrium cases are characterized by the same type of relevant excitations. 
These grow or shrink into the metastable sea depending upon the competition between their surface, which hampers 
cluster growth, and their bulk, which favors it. 

The nucleation and growth of a cluster is controlled in equilibrium, p = 0, by the zero-field interfacial free-energy, or 
surface tension. For p > 0, even though a proper bulk free-energy function does not exist, it seems sensible to assume 
that the properties of the nonequilibrium interface will be determined by a sort of nonequilibrium or effective surface 
tension aQ{T,p) ^. This relies on the assumption that the normalization of the probability measure for interface 
configurations can be interpreted as a nonequilibrium partition function. More formally, let y denote an interface 
microscopic configuration — e.g., the set of heights yk, k = 1 . . . L, that characterize the steps of the interface — 
, and Z~^fi{y;T,p) the associated invariant measure, which for < p < 1 is non-Gibbsian. Then, we define 
ao{T,p) = —{j3L)^^\nZ ([25h. Similar hypotheses have already been tested in relation with other nonequilibrium 
phenomena. For instance, concerning nonequilibrium phase transitions, the distribution of complex zeros of the 
invariant measure normalization factor has recently been shown to obey the Lee- Yang picture in the case of the 
asymmetric simple exclusion process (|20h and for systems with absorbing states ([2ll): see also (|2^ ; l25l). Interestingly, 



In general, the surface tension depends on the orientation angle of the interface with respect to one of the lattice axes, ao{(l>,T,p). 
However, the surface tension entering nucleation theory is the one defined along a primitive lattice vector Q[^[^[3),so</) = henceforth. 
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FIG. 3 (Color) Configurations as the system escapes from tlie metastable state, at different times as indicated, for L = 128, 
p = 0.01, h = -0.25, and, from left to right, T = 0.1, 0.3 and 0.7 in units of To„s. 



the normalization Z is defined in general up to a multiplicative factor (20). This is related to the fact that one has 
to arbitrarily choose a reference interface configuration from which to measure the relative probabilities of all other 
microscopic configurations. However, in the case of our model, the ambiguity in Z can be easily resolved by noting 
that a well-defined equilibrium limit exists, p = Q. In this way, demanding Z to converge to the right equilibrium 
limit for p — Q, the ambiguity disappears. In more general nonequilibrium cases, where no proper equilibrium limit 
can be defined, it is expected that the spurious multiplicative factor won't play any relevant physical role ( 20). 

Computing analytically the exact interfacial invariant measure is beyond our goals, but a solid-on-solid 

approximation ([26l ) will suffice. It thus follows (we refer the reader to Ref. ((25l ) for further details) that ao{T,p) has 
non-monotonous T dependence for any p > and a maximum at a non-trivial value of T which depends on p, see 
Fig. m This peculiar behavior can be understood in simple terms. From low- and high-T expansions of (7o{T,p), one 
obtains that 



— TTT- (Jo + small corrections , 



T. 



(I) 



(9) 



of 



where T^\t,p) is an effective interfacial temperature and crdT) is the exact equilibrium surface tension ([23). At 
high— T, where thermal fluctuations dominate over the nonequilibrium perturbation, T^P is proportional to T — see 
inset to Fig. |4] — , so that tro is proportional to CTo which decreases with T. However, xj^p saturates to a constant as 

T ^ 0, namely, limr-^o T^j\T,p) = 2/ ln[(l - ^)/{p + y/p)] (HI)- This, which is due to non-vanishing nonequilibrium 
fluctuations, implies that ctq is a linearly-increasing function of T at low-T. Consequently, it follows a non-monotonous 
T-dependence of ctq with a maximum which roughly agrees with the crossover observed in T^j^ . 

The properties of the nonequilibrium interface may be obtained from the resulting ao{T,p). For instance, the 
temperature dependence of this function leads to the phase diagram Tdp). In equilibrium, the interface free energy 
goes to zero as T increases toward Tons, and there is no surface tension in the disordered phase for T > Tons (|28l ). 
Therefore, one may identify Tc{p) as the (finite) temperature for which ao{T,p) = 0. We thus obtain the result in Fig. 
[TJ On the other hand, (Tq is linear in T as T — > (see Fig. U), namely cto(T,p) « a{p)T, with (f23) 



a{p) — In 



p + Vp 



(10) 



The condition a{pc) = thus signals the onset of disorder at low temperature. This yields Pc = (\/2 — 1)^ w 0.1716, 
in excellent agreement with previous MC simulations (fill ; [l9| ) . 
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FIG. 4 (Color) The nonequilibrium surface tension as a function of T for, from top to bottom, p = 0, 10"'', 10'^, 10"*, 10"^, 5 x 
10 ,10 , 2 X 10 , 3 X 10 , and 4 x 10 . For any p > 0, the surface tension behaves non-monotonously, contrary to the 
equihbrium case. Inset: The effective interface temperature T^\t,p), as defined in the main text, versus T for the same values 
of p than in the main graph. Notice that T^\t,p > 0) strongly deviates from T in the low-temperature regime. 



C. Effective Free-Energy and Nucleation 

The structural similarities between equilibrium and nonequilibrium excitations responsible for the metastable-stable 
transition allows one to write an effective free energy associated to a cluster of radius R. In formal analogy with the 
equilibrium case (Q), we assume that 

TiR) = "f [2QRao - nR'^2mo\h\] . (11) 

Here, do is the nonequilibrium surface tension in the previous subsection, mo{T,p) is the zero-field spontaneous 
magnetization, which follows from the mean-field scheme in ijlll.Al and fl{T,p) is the cluster form factor, which 
relates the cluster linear length R to its volume, C — flR^. This may be computed from ao via the Wulff construction 
([25I : [29I ). The multiplicative factor in eq. (|lip . 7 ~ 1, is a phenomenological parameter, which might have a very 
weak dependence on T and p, intended to capture small corrections to classical nucleation theory due to correlations 
between different clusters, inhomogeneous magnetization profiles within a cluster, etc (0)- 

The first term in J-{R) measures the cost of the cluster interface, while the second term stands for the bulk gain. 
In this way, small clusters — with a large surface/volume ratio — tend to shrink, while large clusters tend to grow. 
The critical cluster radius, 7?,c, which corresponds to the maximum of J-{R), separates these two tendencies. One has 
that 

We measured this quantity in MC simulations. With this aim, one produces an initial state with a single square cluster 
of down spins and side 2R (the stable phase) in a sea of up spins. This is highly unstable, so that any subcritical initial 
cluster, R < TZc{T,p), will very quickly shrink, while supercritical ones, R > TZc{T^p), will rapidly grow to cover the 
whole system. Since our dynamics is stochastic, we define the probability that a cluster of size R is supercritical, 
Pspc{R). This is measured in practice by simply repeating many times the simulation and counting the number of 
times that the initial cluster grows to cover the system. TZc is then defined as the solution of PspdT^c) = 0.5. As 
illustrated in the inset to Fig. [5l Pspc{R) shows a sharp transition from to 1, which allows for a relatively accurate 
estimate of TZc- ^ The agreement shown in Fig. [5] between these results and our analytical predictions is excellent 
for temperatures well below Tc{p). Interesting enough, TZc{T,p > 0,h) exhibits non-monotonous T-dependence, as 
expected from the anomalous low-T behavior of cto- 



^ Our numerical results on TZc depend in principle on the shape of the initial excitation, which we choose squared. However, we checked 
that using spherical initial excitations (good at high T) one obtains estimates for TZc which are very close to those obtained with squared 
excitations (expected at low T, see Ref. ||25D ). 
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FIG. 5 (Color) TZc vs. T"^ for L = 53, /i = -0.1 and, from top to bottom, p = 0, 0.001, 0.005 and 0.01. Symbols are MC 
results. Lines correspond to the theoretical prediction. The n*'' curve (from bottom to top) has been shifted (4 — n) units in 
the y-axis. Inset: MC results for Pspc{R) as a function of R for L = 53 at T = 0.4 Tons, P ~ 0, and h = —0.1. All the data in 
these graphs correspond to an average over 1000 independent experiments. 
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FIG. 6 (Color) Cluster distribution P(C) for T = 0.2 Tons, L = 53, h ^ -0.1 and, from bottom to top, p = 0.001, 0.005, 0.01, 
0.02, with 7 — 0.815, 0.82, 0.83, 0.85, respectively. Lines are theoretical predictions and points are MC results. The n*'' curve 
(from bottom to top) has been rescaled by a factor 10'"~^' in the y-axis. 



Our ansatz in eq. pT|) also implies that the fraction of stable-phase clusters of volume C = ilR^ in the metastable 
phase follows a Boltzmann distribution 



P(C) = exp ^-2/37 l^aoVnC - mo|/i|CjJ (13) 

Here, the normalization A4 = 28/(1 — m), where Q = Ymj=i ^ 6xp[— /9JF(C)] with C* = riT?.^, is defined so that the 
metastable state has magnetization m{T,p,h) as derived in i jlll.AI The cluster distribution P{C) which is obtained 
in MC simulations for times much shorter than the relaxation time is shown in Fig. [S) This compares very well 
with the analytical prediction, eq. ()13p . For p = 0.001, we observe a non-trivial structure in P{C) which is not 
captured by our continuous description. This is due to the underlying lattice anisotropy, which for low-T and small 
p gives rise to resonances in P{C) for clusters with complete shells, i.e. 2x2,3x2,3x3; see Fig. [6l For larger p 
and/or T, fluctuations smear out these resonances. Notice in Fig. [6]that the nonequilibrium perturbation p enhances 
fluctuations and favors larger clusters. Similar good agreement for P{C) is found in a wide range of values for T, p 
and h. On the other hand, systematic corrections in P(C) show up for microscopic clusters. These appear because 
our main ansatz, eq. pT|) . is based on mesoscopic considerations and therefore is only valid for large enough clusters. 
The observed deviations won't play however a relevant role in our discussion, as far as the critical cluster remains 
large enough, as it is the case here, see Fig. [5] 

The assumption above on T{R) implies a force per unit area which controls the growth of supercritical clusters. It 
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FIG. 7 (Color) Propagation velocity vq vs. for h = -0.1 and p = 0, 0.001, 0.005, 0.01, 0.02, 0.03, 0.04, 0.05 (from 

bottom to top). Curves are linear fits to data. Inset: Slope of linear fits vs. p. The line is the theoretical prediction. 



follows that the propagation velocity of a cluster in the large-size limit will obey the AUen-Cahn expression 13; 0): 

.0 = ^^^, (14) 

where F' is a non-universal constant. We measured vq in MC simulations of an infinitely-large cluster. This was 
implemented in practice by a flat interface separating the stable phase from the metastable one which constantly 
propagates. The results, which are shown in Fig. [71 agree with our predictions. We also observe that vq exhibits 
non-monotonous T-dependence, as expected given the behavior described above for cto. This implies, in particular, 
that, at low T, cooling the system favors the interface propagation or, in other words, that the interface dynamics 
becomes more sluggish as we raise the temperature in the low-T regime. 

As discussed before, ctq increases linearly with T at low temperature, namely, ctq ~ Oi{p)T as T — > with a{p) given 
in eq. pO)). Therefore we expect in this regime that wq ~ [r/ct (p)] {Tons/T) with F = 2r'|/i|rQ~g assuming that 
mo{T 0,p) ^ 0(1). This low-T behavior, linear in T^^, is also confirmed by simulations, see the inset to Fig. [3 
with F « 0.0077. 

The nucleation rate X, i.e., the probability that a critical cluster nucleates per unit time and per unit volume, may 
be written as 

I = A\hfexp[-pJ^inc)], (15) 

where A{p) is a non-universal amplitude and (5 w 3 for random updatings (6). Two main relaxation patterns arise 
depending on the interplay between the relevant length scales in the problem. In fact, one needs to deal with the 
competition between the system size, L, and the mean cluster separation, TZq{T,p, h). The two other important length 
scales, namely the critical radius, TZc, and the thermal correlation length in the metastable phase, fms, do not compete 
with the former two given that they are much smaller in the most interesting regime, ^,ns <C TZc ^ Ti-o, L. Now TZq can 
be calculated together with its associated time scale by requiring that 7?.o = vqIq, which states that a supercritical 
cluster propagating with velocity vq will grow a distance H.^ in a time Iq, and that TZqIoT = 1, which states that the 
probability of nucleating a volume TZq in a time to is (by definition) one. Therefore, TZq = (vo/^Y^^- For TZq ^ L, the 
random nucleation of a single critical cluster is the relevant excitation. In this regime, known as Single-Droplet (SD) 
region, the lifetime of the metastable state is tsd = {L^^)~^ . On the other hand, for TZq ^ L, the metastable-stable 
transition proceeds via the nucleation of many critical clusters, a regime known as Multi-Droplet (MD) region, and 
TMD = [f^z;gl/(31n2)]-i/3 Summing up, 

(rr M / (^'^)"' 7^o»L(SD), 

^ ^ \ (v^l) 7^o«i(MD). ^ ' 

The crossover between these regimes, given by the condition TZq = L, characterizes the dynamic spinodal point, 
namely, 



(17) 
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FIG. 8 (Color) Lifetime r of the metastable state vs. T for the same values of p and conditions than in Fig. [5] The n"* 
curve (from top to bottom) is rescaled by a factor 10"^*^""^^ Solid lines are theoretical predictions for (from top to bottom) 
7 = 1, 0.85, 0.77, 0.65. Amplitudes A{p) are in the range [10"'^, 10"'^]. Dashed lines are predictions for 7 = 1. 
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FIG. 9 (Color) The dynamic spinodal !/idsp|, as defined in eq. gTl, versus T for p = 0, 0.001, 0.005, 0.01, and 0.02 from top to 
bottom, respectively. We explicitly indicate the T ranges corresponding to the SD and MD regimes in the case of the example 
in Fig. m namely, for \h\ = 0.25 and p = 0.01. The symbols x mark the temperatures T/Tons = 0.1, 0.3 and 0.7 for \h\ = 0.25. 
Notice that, unlike at equilibrium, there is here a low-T regime where the MD scenario holds. This is confirmed in Fig. [31 



That is, one should observe the SD regime for \h\ < |ft.DSp|, while the MD regime is expected to dominate otherwise. 

We estimated r from MC simulations by defining it as an average of the first-passage time to zero magnetization. 
Due to the strong local stability that characterizes the metastable state at low temperature, r may be extremely large 
in practice. (For example, this goes up to 10""^ MC step per spin in Fig. [51 Assuming the MC time unit corresponds to 
the inverse of the typical phonon frequency'^, which is of order of 10~^^ s, this corresponds to t ^ 10^^ s!). Therefore 
we needed to use an efficient rejection-free algorithm. Our choice was the s — 1 variant of the MC algorithm with 
absorbing Markov chains, together with the slow-forcing approximation (jsil ). Fig. [5] shows our results for r vs. T 
and different values of p, together with theoretical predictions. Interesting enough, we observe that r increases with 
T at low-T for fixed p > 0. That is, the local stability of nonequilibrium metastable states is enhanced at low-T 
as the strength of the thermal noise (T) increases. This, which is in contrast with the Arrhenius curve observed in 
equilibrium, closely resembles the noise-enhanced stability (NES) phenomenon reported in experiments on unstable 



^ Certain rare-earth magnetic materials can be modeled to a good approximation by the Ising Hamiltonian (30). When we talk here 
about the phonon frequency wo refer to the phonon frequency in one of these real magnetic system with Ising-like properties. 
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systems ((s^). On the other hand, increasing p for fixed T always results in shorter r. This complex phenomenology is 
captured by our simple ansatz eq. Ijlip . which traces back the NES phenomenon to the low-T anomaly in ctq, which 
in turn is a nonlinear cooperative effect due to the interplay between thermal and nonequilibrium fluctuations. That 
is, though any of these two noise sources (T and p) will separately induce disorder, their combined effect results in 
(resonant) stabilization of the nonequilibrium metastable state. 

This interesting effect is also evident in the morphology of the metastable-stable transition. In particular, |/idsp| 
inherits the non-monotonous T-dependence of coj see Fig. [HI resulting in a novel MD regime at low-T not observed 
in equilibrium. For instance, for h = —0.25 andp = 0.01, we expect from Fig. [S]MD morphology for, e.g., T = 0.1 To„s 
and T = 0.7 Tons but SD behavior for T ~ 0.3 Tons, which is confirmed by direct inspection of escape configurations 
such as the ones in Fig. [31 Notice that the SD decay mechanism is a finite-size effect, and that MD will be the 
only relevant mechanism for large enough L. However, the relevance of the SD regime vanishes logarithmically with 
L — see eq. PT|| — so that it could be observable in mesoscopic samples, at least (Q). 



D. Langevin Type Description 

We may get further insight into the above stochastic resonance phenomena and the peculiar low-T behavior by 
rewriting the rate ([2]) as 

p+(l-p)*(/3A7^,) = *(/3efAH0. (18) 
The resulting effective temperature Td = P^^^ then follows as 



Tef(T,p; AH.) = ^ (19) 



AT-Ci+Tln 



For any p > 0, T^f changes from spin to spin and depends on the local order as characterized, for instance, by the 
number rji of neighbors that are in the same state than the central spin (in fact, one can write ATii ~ 4(?7i — 2) + 2sih, 
see Table I): the larger rji, the higher the amount of local order, and the larger the local effective temperature. This 
observation is not restricted to the Glauber function we use for in this paper, see §11 and eq. but holds for 
many other known forms for 5*. For instance, for the Metropolis case, which is analytically simpler for the purpose, 
one obtains for /i — > that 

f -8{\n[p+{l-p)e-^('])~' 77, =4, 
T^:"\t,p)^1 -4(ln[p+(l-p)e"4/3])-i ^^^3, (20) 

[t ??» = 2,1,0, 

such that T^f^ > T^^p > yj^^'^'"^ > Q, the equality holding only in the equilibrium limit p = 0. Therefore, in all 
cases, the strength of fluctuations that affect a spin increases with the local order parameter. This is the fingerprint 
of multiplicative noise, which suggest us to write a Langevin equation with the hope it will capture the physics of 
the problem. Defining the coarse-grained order parameter as the local magnetization field ip(r,t), we propose the 
following Langevin equation 

dij{r, t) 



dt 



DV'^tjjir, t) + V'(r, t) - V'(r, tf + h+ y/ D + fiiP{r,t)^^{r, t), (21) 



where ^(r, t) is a Gaussian white noise with (^(r, t)) — and (^(r, t)£^{r' , t')) — 2S{t — t')S{r — r'), D is the strength of 
the thermal (additive) noise, h is an external field parameter, and /i is the renormalized version of the nonequilibrium 
parameter p. This equation describes the coarse-grained dynamics of a Ising-like system in the broken-symmetry 
phase under an external magnetic field, with a noise term whose amplitude increases as the local order parameter 
ip{r, t) increases, in close analogy with our observations above. Solving the full spatially-dependent problem, eq. (j21|) . 
is a complex endeavor which we will not undertake here. However, as a proof of concept, one can analytically solve a 
simplified model corresponding to the mean-field, 0-dimensional form of the above Langevin equation: 

9-0 



_^^^_^3^h+ ./DTWm , (22) 

where now '0(t) has no spatial dependence. Despite the obvious simplifications, this equation still contains the essential 
competition between thermal {D) and nonequilibrium (/i) fluctuations in a metastable potential that characterizes 
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our system . In particular, eq. (|22p describes the motion of a Brownian particle in an asymmetric bimodal potential, 
namely, V{ip) — — ^V'^ + ~ ^"0 — which is depicted in the inset to Fig. [TOl The particle is subject to fluctuations 
which increase with ^p'^ and whose amplitude remains non-zero as D ^ for any /i > 0. Notice that the noise 
strength ^/D + fMjj^ increases with the order parameter (~ ■0^), in accordance with our observation above concerning 
eqs. ([T ^ - (|20p . For /i = 0, this reduces to the additive theory which describes simple symmetry-breaking systems (j33l ) 
while, for > 0, additive and multiplicative fluctuations compete in eq. (j22p . 

The Stratonovich stochastically-equivalent Fokker-Planck equation of eq. is dtP{ip,t) — —d^[A{il;)P{ilj,t)] + 
ia^[S(V')P(V',0] with A{ip) (1 + - V'^ + ^ and B{ip) = 2{D + fiip^). This yields the steady distribution 



A 



exp 



D~ 



where A is a normalization constant, and 



(23) 



W{ip) = — 




(24) 



The extrema of the effective potential W{^) are the same than those of V^('(/'), namely ipk — 2 cos(0fc)/-\/3, with 



^fc = |[cos ^{—\^^2^K) + 2A:7r], and /c = 0, 1, 2, see inset to Fig. [TOland Fig. [TTJ For < 0, V'o, V'l and ^2 correspond 



to the metastable, stable and unstable extrema, respectively, and the escape time from the metastable minimum is 



t(I?,^,/i) 



27r 



^\V"{^o)V"{^2)\ 



exp 



(25) 



where V" = d^V{x). Identifying D with the temperature parameter in our microscopic model, this precisely repro- 
duces the thermal NFS phenomenon for r as described above; compare Figs. [TO] and [5] Therefore, the resonant, 
temperature-enhanced stabilization of nonequilibrium metastable states may be regarded as a consequence of com- 
petition between additive and multiplicative noise. 



E. The Limit of Metastability 

When the magnetic field \h\ is increased, the strength of the metastable state decreases. That is, the local minimum 
in the free-energy functional associated to the metastable phase becomes less pronounced. Upon further increasing \h\, 



As a result of the homogeneity hypothesis associated to the mean-field approximation, eq. H22| l has no spatial structure and therefore 
is unable to describe the compact excitations responsible of the metastable state decay in the microscopic model. Despite this obvious 
difficulty, eq. I|22| l still contains the relevant competition between additive and multiplicative fluctuations which gives rise to resonant 
behavior in this model. 
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FIG. 11 (Color) Effective potential W{4>, D, fi) for fi — 4 x W ^ and h — —0.1. Note how the energy barrier between the two 
local minima changes as a function of D for this fj, > 0. 



such local minimum eventually disappears. At this point, the metastable state becomes unstable, which means that 
its relaxation toward the final stable configuration is no longer hampered by any free-energy barrier. The magnetic 
field strength h*{T,p) at which this metastable-unstable transition occurs is known as the spinodal {H) or, sometimes 
(see Appendix A) pseudospinodal or intrinsic-coercive field. 

As \h\ is increased, the metastable phase remains homogeneous on a coarse-grained length scale, for times shorter 
than T, so that the mean-field approach in section UlI.AI is suitable within this context. The locally-stable steady 
magnetization, as computed from the stationary solution of eqs. (|H) and ([5]), exhibits two branches as a function of 
h, see the inset in the left graph of Fig. [121 This hysteresis loop reveals that, at the mean-field level, metastability 
disappears abruptly for any \h\ > h* . In order to evaluate h*(T,p), one may study how the metastable state responds 
to small perturbations of the applied field. If {x^^,z^^) is a locally-stable stationary state for T, p and Hq, with 
magnetization opposed to Hq, and we perturb h = ho + Sh, the new locally-stable stationary solution is modified 
according to x^^ = + and z^^ — + ■ One obtains to first order that 



dhG2d.Gi - dHGid,G2 



d^Gid,G2 - d^G2d,Gi 



xSh, (26) 





and a similar expression for e^, where the quantity in brackets is evaluated at {x^^l',z^l') for given values of the 
parameters T, p, and ho, and the functions Gi^2ix,z) are defined in eqs. ([4]) and ([5]). This (linear) response diverges 
for 



[d:,Gid,G2 - a,G25.Gi]o = , (27) 

which corresponds to a discontinuity in the metastable magnetization as a function of h. For fixed T and p, the 
field for which eq. (|27p holds is h* . The left graph of Fig. [12] shows the mean-field result for h* (T, p) as obtained 
numerically from eq. (|27p . In particular, for p — 0, we recover the standard equilibrium mean-field spinodal curve: 
Converging to 2 as T — > 0, linearly decreasing with temperature for small T, and vanishing as (Tbothc — T)^/^ at the 
mean-field equilibrium critical point. For p > 0, an instability separates two different low-T behaviors of h*{T,p) 
depending on the amplitude of the nonequilibrium fiuctuations. For small enough values of p, namely, p € [0, 0.031], 
which includes the equilibrium case, the field h*{T,p) monotonously grows and extrapolates to 2 as T ^ 0. For larger 
p, up to p = 5/32, however, h*(T,p) — > as T — > 0, with a maximum at some intermediate value of T. 

The value p — ttc ~ 0.0315 separates these two types of asymptotic behavior. The p < tTc regime can be understood 
on simple grounds. In this case, h*{T,p) increases as T drops, meaning that a stronger field is needed to destroy 
metastability. In a metastable state, the tendency to maintain spin order prevails over the tendency of the individual 
spins to follow the external field. Since both T and p induce disorder, one should expect that decreasing T and/or p 
a stronger field will be needed to destroy the metastable state, as confirmed in Fig. [T^lfor p < tTc- 
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FIG. 12 (Color) Left panel: Mean-field result for h*{T,p) as a function of T for, from top to bottom, p = 0, 0.01, 0.02, 
0.03, 0.031, 0.032, 0.035, 0.04, 0.05 and 0.1. The qualitative change of behavior in the low temperature region occurs for 
p G (0.031,0.032). We also indicate for \h\ = 0.25 the temperatures Ti < T2 which, as discussed in the main text, comprise 
the metastable region for p — 0.05. The inset shows the two locally-stable steady magnetization branches as a function of h 
for T — 0.77^(0) and p = 0.005. The solid (dashed) line represents stable (metastable) states. The dot-dashed line signals the 
discontinuous transition, at /i = h*{T,p), where metastable states disappear. Right panel: Monte Carlo results for h*{T,p) 
as a function of T for L = 53 and, from top to bottom, p = 0, 0.01, 0.03, 0.0305, 0.0320, 0.0350, 0.04 and 0.05. The change 
of asymptotic behavior in the low-T limit happens for p £ (0.03,0.0305). The inset shows the probability of occurrence of 
metastability, as defined in the main text, as a function of ft < for L = 53, T = O.JTons and p = 0. Data here correspond to 
an average over 500 independent demagnetization runs for each value of h. 

The p > TTc regime is more intriguing. Consider, for instance, p — 0.05 > tTc and \h\ = 0.25. According to the left 
graph in Fig. [T^ we can define two temperatures such that metastability only occurs for T S [Ti, T2]. The fact that h* 
goes to zero as T ^ in this regime means that, at low temperature, the nonequilibrium fluctuations are strong enough 
to destroy on their own the metastable phase. Based on the above argument, one would expect that increasing T 
should avoid metastability. However, there is metastability for intermediate temperatures, T g [Ti, T2]. This reentrant 
behavior of the spinodal field suggests once more a resonance between thermal and nonequilibrium fluctuations: despite 
the intrinsic disordering effect of both noise sources, T and p, they cooperate to produce metastability, for which local 
order prevails. We measured h* in MC simulations, confirming this; see the right graph in Fig. 1121 and appendix A. 
In particular, we find the instability for p — tTc ~ 0.03025 in simulations, which is very close to the mean-field result. 

There is an important difference between the two resonant phenomena above, i.e., the noise-enhanced stabilization 
of nonequilibrium metastable states and the reentrant behavior of the spinodal field. That is, while the former occurs 
for any p > 0, the latter only occurs for p > tTc- In order to understand this fact, we studied the limit of metastability 
within the cluster nucleation scheme of section [ill. CI In this case, h*(T,p) may be estimated from the condition 
that the critical cluster contains one single spin, 2TZc(T^p, h*) = 1. This is equivalent to require that any microscopic 
fiuctuation can trigger the transition toward the stable state, which signals the onset of instability. The result, namely, 
h*(T,p) — ctq/too according to eq. (|12p . must be taken with some care, since it assumes that the surface tension cto, a 
macroscopic concept, is relevant for microscopic clusters. In fact, this corresponds to a mean-field spinodal Q), which 
exhibits reentrant behavior for any p > 0. Notice that the nonequilibrium perturbation in our model is stochastic in 
origin, producing random spin flips with probability p. For a macroscopically large interface the probability that these 
stochastic perturbations do not affect its structure and dynamics is negligible, so the nonequilibrium perturbations 
are always present, inducing the observed low-T anomaly in ctq, see ijlll.BI On the other hand, microscopic clusters 
(as those involved in the computation of h* ) have a non-negligible probability of not being affected by the stochastic 
perturbation p during their nucleation and evolution for small enough p. In this way, these small clusters unaffected 
by p would behave effectively as equilibrium clusters. This would renormalize the (microscopic) critical cluster size 
to a new, larger effective value, which thus would modify our estimation of the spinodal fleld h*. Hence one expects 
h* to behave as in equilibrium for small enough values of p, as is in fact observed, while for large values of p (say 
p > TTc) even small clusters are affected by the nonequilibrium perturbation, involving reentrant behavior for h*. The 
challenge remains to put this qualitative argument, which seems plausible, into more quantitative grounds. 
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FIG. 13 (Color) Left: Schematic representation of the model particle with open circular boundary conditions. Spins at the 
border do not have nearest neighbor spins outside the circle. Right: Typical evolution in which the magnetization is observed 
to decay by jumps to the final stable state. This is for a single particle of radius r = 30 (~ 10^ spins) at low temperature, 
and for small values of h and p (see the main text). The time axis shows t — tq in MCSS (Monte Carlo steps per site) with 
To — lO'^^MCSS; this is of order of the duration of the initial metastable state. The inset shows a significant detail of the 
relaxation. 

IV. SMALL PARTICLES 

The above concerns bulk metastability, i.e., a property of macroscopic systems or models in the thermodynamic 
limit. However, the case of finite systems, in which surface effects are relevant, is of great interest. This occurs for 
magnetic particles in many applications, e.g., magnetic recording media where one often needs to deal with particle 
sizes ranging from mesoscopic to atomic levels, namely, clusters of 10^ to 10^ spins. We therefore also studied our 
model with open circular boundary conditions. The system is now defined on a square lattice, where we inscribe a 
circle of radius r, and assume that the bonds to any spin outside the circle are broken; see the left graph in Fig. [T5] 

The effects of free borders on the metastable-stable transition have already been studied in equilibrium systems 

(0; 113; M)- 

In this case, the system evolves to the stable state through the heterogeneous nucleation of one or 
several critical clusters, which always develop at the border, where it is energetically favorable. Apart from this, 
the properties of the metastable-stable transition do not change qualitatively as compared to the case of periodic 
boundaries. In our nonequilibrium system we also observe this kind of heterogeneous nucleation. However, in closer 
inspection, one notices that the randomness of the decay process is importantly enhanced, which results in some 
unexpected phenomenon. 

The main observation is illustrated in the right graph of Fig. 1131 That is, when observed on the appropriate time 
scale, namely, after each MC step (per spin), the time relaxation of the magnetization, m (t) , from the metastable 
state occurs by strictly monotonic changes. There is thus a sequence of well-defined abrupt jumps of m{t), which 
we will name avalanches in the following, that resembles familiar relaxation processes in many complex systems in 
nature. We checked that this is a general feature of the model relaxation at low T. However, a too rapid relaxation 
and/or domains too diluted and interfaces too diffuse will tend to obscure the observation. Consequently, it turned 
out preferable to deal with small values of temperature, in order to have compact and well-defined clusters, and with 
values of h that do not excessively accelerate the evolution. In addition, the parameter p can take a wide range of 
values, provided that its effects are comparable to those of temperature. A perfect compromise is T = 0.11 Tons, 
p = 10~^, and h = —0.1, and most data below concern these parameter values. 

A. Avalanche Statistics 

To be precise, consider an avalanche that begins at time ta, when the system magnetization is m{ta), and finishes 
at tb- We define the avalanche duration as At = |ia — ^bl and its size as = \m(ta) — m{ti,)\, and we are interested on 
the associated respective distributions, P (At) and P (Am). It turns out to be most important for the reported result 
to remove from the data some trivial extrinsic noise (|35r). namely, small thermal events of typical size (see appendix 




FIG. 14 (Color) Left: Log-log plot of the size distribution P (Am) of avalanches (i.e., the data after subtracting "small events" 
as defined in the main text) for ensembles of independent model particles of radius (from bottom to top) r = 30, 42, 60, 84 and 
120, respectively. This corresponds to the O symbols (green). Plots of the corresponding duration distribution -P(Ai) vs. cA^ 
for each r are also shown (x, red) superimposed, with c ~ 0.5 and 7 ~ 1.52 (see the main text). For visual convenience, the 
curves are shifted vertically by 4" with n = to 4 from bottom to top. Right: Log-log plot of the duration distribution P (At) 
for the same ensembles. For visual convenience, these curves are shifted vertically by 2" with n = to 4 from bottom to top. 
The inset shows a log-log plot of the size (top) and duration (bottom) cutoffs A* vs. r. Lines are power-law fits. Running 
averages have been performed in all cases for clarity purposes. 
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(28) 



These events correspond to the short-length fluctuations that are evident by direct inspection in the inset of the right 
graph in figure [T31 These also correspond to avalanches that originate at flat interfaces which, at low enough T, have 
a significant probability to form, as discussed below. 

The distribution P (Am) that results after deducting these small events is illustrated in the left graph of Fig. [TH 
This is well described by 



P(A„,)~A,-/W, 

with a size-dependent exponent given by 

T (r) = Too + air"^ , 

where Too — 1.71 (4). The left graph in Fig. [14] depicts the corresponding duration distributions. This follows 



— a(r) 



(29) 
(30) 

(31) 



with 



a (r) ^ aoo +a2r ^ , (32) 
where aoc = 2.25 (3). In both cases, size and duration, the apparent power law ends with an exponential tail, i.e., 

P(A) -exp(-A/A*) . (33) 

The cutoffs that we observe follow A* ~ with Pm ~ 2.32 (6) and ~ 1.53 (3), respectively (see the inset in the 
left graph of Fig. [T4)) . 

In practice, we observed this behavior for several combinations of values of the parameters T, p and h. The 
requirement for an easy observation is that, as indicated above, the configurations need to be sufficiently compact and 
that none of the underlying processes should predominate so that it obscures the other. That is, the above behavior 



17 




FIG. 15 (Color) Typical time evolution of a circular model particle of radius r = 120 as it decays from the metastable state. 
Time increases from left to right and from top to bottom. The metastable (stable) phase is the blue (red) region, and we 
represent the specific avalanches at each time as yellow regions. As discussed in the main text, all the large avalanches (the 
only ones that are visible here to the naked eye; the small others are indicated by arrows) occur at the curved parts of the 
interface. This simulation is for T — 0.11 Tons, p ~ 10~^ and h = —0.1. 

is the consequence of a competition between metastability, thcrmalization and the nonequilibrium perturbation. As a 
matter of fact, as it will be discussed below, observing power laws requires both free borders and the nonequilibrium 
condition. That is, the distributions P (A) look approximately exponential if the system has periodic borders and/or 
one sets p = in (21). It is also obviously convenient, in order to have good statistics, a choice of the parameter values 
that selects a region of the phase diagram in which evolution proceeds by as many jumps as possible. In any case, 
however, the behavior that we describe here is robust and characterizes our model, and no fine-tuning of parameters 
is needed. 



B. Many Scales 

We next show that the observed power laws are in the model a superposition of many individual processes, each 
with a well-defined scale, and that these may be traced back to specific demagnetization events. Therefore, there is 
no need here to invoke for any critical point, which is the familiar origin for scale invariance in equilibrium. 

Let us consider a relatively-Zarge circular model particle. Inspection of configurations then clearly reveals that, at 
sufficiently low T, the interesting events always occur at the interface between the metastable and stable phases. One 
observes curved domain walls due to the faster growth of the stable phase near the concave open borders. This is 
in accordance with a previous conclusion that, in equilibrium, the critical cluster always nucleates at the free border 
(0; ITsI ). Any curvature costs energy, so that the large avalanches tend to occur at curved domain walls, which will 
tend to produce flatter interfaces; see Fig. 1151 This is confirmed when one monitors the mean avalanche size (Am) 
and the wall curvature (k) as a function of magnetization m {t) . We define the curvature k as the number of ascending 
steps at the interface. That is, the number of up spins which have two ups and two downs at the respective sides 
along the interface. (This definition requires well-defined compact clusters, which are the rule at low T.) The left 
graph in Fig. illustrates our results. After averaging over many runs, definite correlations show up. That is, as one 
could perhaps have imagined, the event typical size is determined by the curvature of the interface region at which 
the avalanche occurs. 

This is also confirmed by monitoring P (A^ | k), the conditional probability that an avalanche of size A™ develops 
at a wall of curvature k. We studied this in detail by simulating an interface of constant curvature that evolves 
according to the transition rate ((H); see appendix C. The right graph of Fig. [TOl shows that P {Am \ n) has two 
regions for given k. The first one corresponds to the small events of typical size ((28)) mentioned above. There 
is also a region which, contrary to the situation in Fig. I14[ exhibits (stretched-) exponential behavior, namely, 
P (Am I k) ~ exp [- {Am/A„,y] with T] « 0.89. That is, a region of the interface with curvature n tends to induce 
avalanches of typical size Am {i^) • 

This fact turns out most relevant because, due to competition between the randomness induced by free borders 
and the one induced by p in the transition rate ([2]) , the interface tends to exhibit a broad range of curvatures with 
time, as illustrated in Fig. [15] and the left graph in Fig. More specifically, relaxation proceeds via a series of 
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FIG. 16 (Color) Left: Semilogarithmic plot of the mean size (Am) (solid, green line) and mean curvature (k) (dotted, red line) 
as a function of magnetization. Notice the non-trivial structure uncovering a high degree of correlation between the avalanche 
mean size and the average curvature of the wall at which it originates. Right: Semilogarithmic plot of P (A^ | k), the size 
distribution for avalanches developing at a wall of constant curvature, k, where k increases from bottom to top (see appendix 
C). Here, Nm = ^NAm- For visual convenience, the curves are shifted vertically by 10" with n = to 4 from bottom to 
top. This corresponds to an average over 3500 independent runs, and running averages have also been performed for clarity 
purposes. 



different configurations, each characterized by a typical curvature of the interface and by the consequent typical form 
of the critical cluster which induces the avalanche. Therefore, what one really observes when averaging over time 
is a random combination of many different avalanches, each with its typical well-defined (gap-separated) size and 
duration, which results in an effective distribution. The fact that this combination depicts several decades (more the 
larger the system is) of power-law behavior can be understood on simple grounds. 

Let Q (A) the probability of A, and P {x \ A) = Aexp {—Ax) the probability of an event of size x given A. Assume 
that A can take a finite number of equally spaced values Aj^, /c = 0, 1, 2, . . . , n, in the interval [^mim ^max]j so that 
Ak — ^min + k5 with 6 — (Amax ^ ^min) /n (alternatively, one may assume randomly distributed ^fes), and that all 
of them have the same probability, Q (A) — const. One obtains that 
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(34) 



The fact that even such a simple, uncorrelated ansatz describes qualitatively the data is illustrated in Fig. [T71 That 
is, the superposition of a large but finite number of exponential distributions, each with a typical scale, yields an 
effective global distribution which is consistent with apparent scale invariance. This distribution extends in practice 
up to a cutoff, exp {—xAjnin) , which corresponds to the slowest exponential relaxation. Notice that eq. ([M]) predicts 
a size-independent exponent, T(r) = Too = 2, that differs somewhat from the observed asymptotic one, Tqo = 1.71(4). 
This reveals that a more complete explanation than ([34]) would be suitable. This would require taking into account, 
for example, dynamic correlations such as the ones revealed in the left graph of Fig. 1161 The above, however, is 
already semi-quantitative and there is no evidence that a more involved computation would modify qualitatively this 
picture. 

Consider next P{At \ A^), i.e., the probability that the avalanche of size lasts a time Aj. As expected from 
the above, this exhibits well-defined peaks corresponding to large correlations, i.e., avalanches of a given size have a 
preferred duration and vice versa. Assuming Am ~ A^, we obtain 7 = Pm/Pt = 1.52(5). Using this relation, one may 
obtain the duration distribution by combining eq. (p4)) with P(Am)dAm = P(At)dA(. A comparison of the resulting 
curve with data in Fig. [T7] leads to 7 ~ 1.52, in perfect agreement with the value obtained from the cutoff exponents 
p. More generally, a scaling plot of P{At) vs. cA] , with c some proportionality constant, must collapse onto the 
corresponding curve P(A„i) for each particle radius r. This is confirmed in the left graph of Fig. [T3]for 7 ~ 1.52, 
further supporting our description of scales superposition. Let us assume for a moment that, as it would occur in the 
presence of a critical condition, both P{Am) and P{At) are true (as opposed to apparent) power-law distributions. 
It then follows the scaling relation (a — 1) = 7(t — 1) which, together with our values above for a and r, implies that 
7 ~ 1.76, which contradicts the value 7 ~ 1.52 that follows by two methods. We believe this misfit simply confirms 
our point that none of the distributions P{A) in this section exhibits true scaling behavior. 
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FIG. 17 (Color) Left: Solid lines are the predictions, equation (|34ll for n = 200, Amin ~ 0.007 and Amax ~ 1. The symbols 
stand for the avalanches duration (lower curve) and size (upper curve) when r = 60, i.e., two of the data sets in figures 2 and 
3. In this particular case, the finite-size exponent is r(r = 60) — 2.06(2). Right: The theoretical avalanche size distribution on 
the left panel, and the many exponentials which compose it. 



C. Some Additional Comments 

We conclude that the model particle relaxes via events, each with a well-defined scale, but many of them randomly 
combine producing (apparent) power-law behavior. This is a consequence in the model of two main features, namely, 
free borders and the nonequilibrium perturbation. (For p = and/or periodic boundary conditions, one just observes 
a well-defined mean.) There is no indication of chaos, e.g., sensitivity to the initial condition. We also discard that 
the observed power laws have any relation in the model with criticality in the familiar, equilibrium sense. In fact, 
we do not observe any relevant correlation other than the dynamic ones that we described in detail. The question, 
given that our model is purposely oversimplified, is whether this picture applies to scale-free fluctuations in natural 
phenomena. Demonstrating this, i.e., analysis of separate elementary events in actual cases, is difficult. However, we 
argue next that there are some indications that this may be the case in some occasions. 

We first mention that our separation between small events — as described in eq. (j28[) and appendix B — and (large) 
avalanches is also supported by experiments (fs^; [sgI ; IstI ). On the other hand, it is remarkable that the statistical 
properties of the resulting distributions are indistinguishable in practice from actual data (38). For instance, size 
corrections similar to the ones in eqs. pop and (|32p for r and a, respectively, have been reported in avalanche 
experiments on rice piles |39l), and our asymptotic values are very close to the ones reported in magnetic experiments, 
e.g. . Too = 1.77(9) , ttoo = 2.22(8) and 7 = 1.51(1) in Ref. (|35r ) for (quasi-two dimensional) ribbons; see also Ref. 



(37). Moreover, our cutoff values in eq. (j33p follow the precise trend observed, for instance, in magnetic materials 
(40; [41I ). Our non-critical scenario is also consistent with the fact that one variously observes either power-law or 
exponential distributions, or a mixture of both, in different but closely related experiments and even in different 
regimes of the same experiment; see, for instance, Ref. (j42( ). 

It is also remarkable that our picture does not imply universality; instead, microscopic details matter, as it is the 
case in experiments. That is, riniversality is suspicious in a context in which critical exponents seem to vary with 
conditions and materials (fssl : l4ll ) , and applications sometimes rely on sensibility to the sample microstructural details 
(|43l:[43). Our picture also explains other features of natural signals such as, for instance, reproducibility. This refers to 
the fact that experiments reveal that the avalanches, i.e., excluding the small events, tend to occur at the same stages 
of evolution (fisl : [i^: IZtI) . As shown in Fig. [181 our system exhibits a high degree of reproducibility due to the strong 
tendency of the critical cluster to nucleate at the border. Simulations indicate that this occurs more markedly the 
lower h and T are, i.e., when the system is more efficient in selecting the most (energetically) favorable configuration. 

Interesting enough, if our picture has a broad applicability, there would be more hope to the goal of predicting large 
events. That is, the assumption of some underlying criticality naturally implies that catastrophes, though relatively 
rare, occur in an strongly correlated bulk and, consequently, have the same cause as the smaller avalanches ('is'; "4^; . 
Instead, in our picture, events are characterized by their size, and each size follows from some specific microscopic 
configuration. The configurations that, under appropriate conditions, may originate large events qualitatively differ 
from the ones corresponding to smaller events. In summary, there is some specific cause for each event which depends 
on its size. 

In spite of the extreme simplicity of the model, one should perhaps recognize that it is likely that its basic features 
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FIG. 18 (Color) The avalanche sizes as a function of magnetization for six typical independent runs. This is to illustrate that 
large avalanches tend to occur, approximately, at the same value of magnetization. This occurs here for the large avalanche 
around m — 0.3 in four of these runs and for the sequence of large avalanches in all the runs for m G [—0.7, —1], for instance. 



are present in some actual systems as well. In particular, one may imagine that the model randomness mimics the 
causes for the random interface rearrangements in magnets and for the slip complexity in earthquakes, for instance. 
Furthermore, in these two cases, as in the model, time evolution is through a series of successive short-lived states 
that constantly halt the system relaxation. It is precisely this what causes a constant change of scale which results in 
the (apparent) scale invariance of the model. 

It would be interesting to study the possible occurrence of these short-lived halt states in natural signals. These 
halts are associated in the model particle with flat interface regions. That is, once the initial metastability breaks down, 
the particle becomes inhomogeneous, and flat interfaces have a significant probability to form after each avalanche 
(which aims at minimizing interfacial energy). As this is the most stable configuration against small perturbations, the 
system remains some time with constant magnetization m (tb) ■ This may be described as an entropic metastability. 
There is no energy barrier but an unstable situation such that a given microscopic random event suffices to initiate 
the next avalanche. 

Finally, we mention that there are many other possible explanations for scale-invariant noise based on non-critical 
mechanisms; see, for example, Refs. (fist Isil: [52l: ISa: Is^: Issi: IsgI) . To our knowledge, however, this is the most general 
one so far reported which has a physics content. Other proposal are often restricted to specific situations and, in some 
cases, they may also be interpreted at the light of a superposition of many different typical scales. A similar origin for 
electronic 1/f noises was suggested in the past (see, for instance, Refs. (H^ [13; \EEi, [H; H^)), though this is perhaps 
the first time in which an explicit relation is drawn between elementary events (avalanches) and microscopic physical 
processes. We also mention that recent observations in a sociological context may be cast into the behavior of our 
"particle". That is, according to certain definitions, there is a relation between fame and merit (|6ll ). It has recently 
been estimated that the probability distribution for a certain level of fame falls off exponentially for a respected 
scientist, while it seems to decay as a power law in other social groups (j63 ) This may reflect that fame in science 
is mainly due to a homogeneous, rather well-informed community, while many different communities, each with a 
different perception of merit contribute to the fame of, for instance, a movies actor. The combination of distributions 
with different means then produces the observed power law. 



V. CONCLUSION 

This paper deals with metastability in a noncquilibrium environment. Wc studied a two-dimensional Ising fcrro- 
magnet subject to competing dynamics, as if a completely random processes were constantly perturbing a tendency to 
thermalization at temperature T. The perturbation impedes equilibrium, even if it is extremely weak, and a noncqui- 
librium steady state sets in asymptotically. This is a convenient background to investigate various questions, and it 
could mimic certain noncquilibrium situations in nature. 

In the light of the observed qualitative behavior, we argue that, at a mesoscopic level, concerning both steady 
and time-dependent properties of clusters in metastable conditions, it is sensible to assume that a noncquilibrium 
free-energy function exists which is formally similar to the equilibrium concept. In particular, we assume that the 
metastable phase results from a contest between surface and bulk terms of this function, and that a noncquilibrium 
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surface tension captures the (strongly nonlinear) interplay between thermal and nonequilibrium noises. We also show 
that a Langevin equation with additive and multiplicative noises captures the essential physics. 

Our theoretical approach predicts resonant coupling among the thermal and nonequilibrium noises which results 
in novel phenomenology. For instance, there is noise-enhanced propagation of domain walls and stabilization of the 
metastable state at low T, and reentrant behavior of the spinodal field under strong nonequilibrium conditions. These 
cooperative phenomena are perfectly confirmed by Monte Carlo simulations. 

We also explored the relaxation from a metastable state in the presence of open boundaries. This case also shows 
intriguing behavior at low T. That is, the decay occurs through a series of many metastable-like states that repeatedly 
halt the dynamics, which resembles the relaxation by avalanches reported for many complex systems, e.g., interface 
rearrangements in magnets and fault slips causing earthquakes. We conclude that scale-free avalanche distributions 
are not to be associated to any critical condition in the model but are simply a superposition of many different classes 
of events, each with a well-defined scale. In the model this is determined by the curvature of the domain wall at 
which the avalanche originates. This may be a property of any complex relaxation phenomena characterized by a 
multiplicity of short-lived, metastable-like states. 
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Appendix A: The Spinodal Field in MC Simulations 

We briefly describe here a method to estimate the spinodal field h* {T, p) which separates metastable from unstable 
states. To be precise, h* is a psewdo-spinodal given that the transition in the actual system (with fluctuations) is not 
sharp but a progressive crossover (jsj), as illustrated in the inset to right graph in Fig. [T^l Estimating h* allows in 
the main text a confirmation of the corresponding mean-field prediction. 

Our method consists in following the system path in phase space as it evolves from the initial state with all spins up 
toward the final stable state. With this aim, we assign to each intermediate state along this path, say s^, a measure 
A(sfe) of the net tendency of the system to evolve toward the stable phase from that configuration. Using the concept 
of spin class described in ^III.Al it is straightforward to see that A(sfe) = G{sk) — S{sk), where 

4 4 

g{sk) = ^ zy(+,n;Sfc)w(+,n); S{sk) = ^ iy{~,n;Sk) uj{-,n) 

n=0 n=0 

are the growth and shrink rates of the stable phase in state s^, respectively (recall that h < 0). Here, v{s,n;Sk) is 
the fraction of spins in class (s,n) for configuration s^, and w(s,n) is the spin-class transition rate. 

Metastability is hampered by free-energy barriers, while unstable states evolve without any impediment. Therefore, 
we may divide relaxation paths in phase space in two sets: Metastable paths, for which at least one configuration exists 
such that A(sfe) < 0, and unstable paths, where A(sfe) > V/c. Given the stochasticity of the dynamics, one needs to be 
concerned with the probability of occurrence of metastability, defined as Il{T,p, h) — rimot/A'oxpj where n-aict{T,p, h) 
is the number of experiments out of the total A'cxp in which the relaxation path in phase space belongs to the class 
of metastable paths. The limit of metastability is defined in this scheme as the field for which Yi[T,p,h*) = 0.5. 
This is illustrated in the right graph of Fig. [12] for a system size L — 53. (We studied finite-size corrections to the 
spinodal field by simulating larger systems; however, these corrections may be neglected for all practical purposes.) 
This confirms the reentrant behavior of h* for low temperature under strong nonequilibrium conditions. 

Appendix B: Small Events 

We now determine the avalanche statistics for a flat interface under a magnetic field at low T. Let us assume a 
macroscopic, perfectly flat interface between the metastable and stable phases which, for ft, < 0, correspond to the up 
and down phases, respectively. 

Flat interfaces are very likely at low T, since they minimize surface tension, and they will tend to invade the 
metastable phase with time for any ft, < 0. 

We are interested in the probability per unit time for a change in the state of a spin in class (s, n), W{s, n), where 
s = ±1 and n G [0, 4] stands for the spin number of up NNs; see i jlll.AI One has that W{s, n) = v{s, n) uj{s, n), where 
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FIG. 19 Schematic representation of a flat interface in which a nucleation event originates and grows as discussed in the main 
text. Notice in (c) that the only relevant spins during this kind of growth are the two spins in class ( — 1,2) and the two spins 
in class (+1, 2) at the respective ends of the cluster. 




FIG. 20 Semilog plot of the avalanche size distributions P{Nm) for avalanches in the field direction (Q) and against the field 
(□), as obtained from Monte Carlo simulations of a flat domain wall of size L = 53 for T = 0.11 Tons, P = 10"® and h = -0.1. 
The slope of the lines correspond to our theoretical predictions in eqs. p5)l and p6)l. 



v{s^n) is the fraction of spins of class (s,n) in the current configuration and Lo{s^n) is the corresponding transition 
rate; see eq. ([2|) and TablelH For the given initial condition, the only possible classes are (+1, 4), (+1, 3), (—1, 1) and 
(—1,0). For < < 2 J and low T, Lo{s,n) is very small for these four classes, e.g., between 10~^ and 1.24 x 10~^ 
for T = 0.11 Tons, P = 10"*^ and h = -0.1. 

Let us assume that a fluctuation nucleates as in Fig. 119b ,. i.e., one spin has been flipped. This spin belongs to class 
(— 1, 3), and its two NN spins in the direction of the interface belong to class (+1, 2). It is straightforward to show that 
the probability that this fluctuation either grows along the interface or shrinks is much larger than the probability of 
any other event in the surrounding bulk. 

Given that there are many more spins in the bulk than in a flat interface, there will be many bulk fluctuations 
before anything happens at the interface. However, whenever an interface fluctuation occurs, an avalanche will quickly 
develop along the wall before anything disturbs the surrounding bulk. 

Therefore we can safely assume that once the interfacial perturbation has appeared, the system dynamics can be 
reduced to the growth and shrinkage dynamics of the interfacial perturbation. Under this assumption, the most 
probable process to be observed consists in the growth of the interfacial fluctuation via the flipping of the lateral spins 
in class (+1,2) at the extremes, until one of the two spins in class (—1,2) which bound the interfacial fluctuation 
flips, halting the avalanche (see Fig. [T9lc). 

We now assume that the system is in the state depicted in Fig. 1191 b. with two up interfacial spins flipped, and we 
want to compute the probability of finding a lateral avalanche of size Nm, i.e. involving Nm spins. The restricted 
dynamics we assume only involves four different spins [two in class (+1,2) and two in class (—1,2)]. Within our 
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FIG. 21 Sketch of the dynamic boundary conditions in the ^-direction. Here we represent the first and last columns of the 
system, and the local, dynamic magnetic field that each spin in these columns feels as the system evolves in time. The effect 
of these dynamic boundary conditions is equivalent to eliminate the interaction of each border spin with its up neighbor, i.e. 
concave open boundary conditions. 



constrained dynamics, the probabilities for the fluctuation to grow or stop are 
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In a similar way we can compute the typical size A^-* of avalanches against the magnetic field, i.e. avalanches 



involving spins in the stable (down) phase. The result is 

A(+) = 
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-20\h\ 



(36) 



As expected, A^^ < aL ^ Fig. shows the size probability distributions P{Nm) for avalanches toward and against 
the field, as obtained in Monte Carlo simulations of a flat domain wall of size L — 53, for for T = 0.11 Tons, P = 10^^ 
and h — —0.1. A comparison with our predictions (|35p and p6p is also shown, with excellent results. 

The small avalanches here described, to be associated with the presence of flat domain walls as the system demag- 
netizes from the initial metastable state, appear in our magnetic nanoparticle together with more structured events 
of large size. As shown in section HVl these structured events exhibit scale-invariant properties which stem from the 
interplay between the nonequilibrium perturbation and the free borders of the magnetic particle. A good measurement 
of these large-scale avalanches involves filtering the above trivial noise, also known as extrinsic noise dsF 



Appendix C: Avalanches from Walls of Constant Curvature 

We have shown in ijlV.BI that large avalanches originate due to domain wall curvature ^. Moreover, domain wall 
curvature appear in the system due to concave open borders, since spins near the concave border flip faster than bulk 
spins. 



In this context, we define the amount of curvature of a domain wall as the number of kinks in the interface, e.g., the number of up 
spins flanked, respectively, by two ups and by two downs at the sides along the interface). This definition requires well-defined compact 
clusters, as for low temperature. 
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FIG. 22 Snapshots of an interface that evolves according to transition rate ^ subject to concave open boundaries as described 
in the text. The stable (metastable) phase corresponds to the green (orange) regions. The yellow regions correspond to the 
avalanches in light yellow. From top to bottom, Lx — 20, 40, and 100. 

The MC simulations reported in §IV.BI show that the size of the avalanche is strongly correlated with the curvature 
of the interface region at which the event originates. Here we deep on this correlation, namely, the interest is on the 
probability that an avalanche of size Nm originates at a wall of some specific curvature. With this aim, we modified 
our basic system so that it shows an interface with constant, except for small fluctuations, non-zero curvature. 

Consider a semi-infinite, x oo lattice with concave open boundary conditions in the x direction. This is done 
in practice starting with a, L^x Ly lattice and fixing the spins at the top row to +1, while spins in the lowest row 
are fixed to —1. On the other hand, boundary conditions in the x direction are dynamic: the lattice is also open in 
the X direction, although each spin in the first and last column suffers an additional dynamic magnetic field, equal at 
any time to the negative value of its up neighbor. This is sketched in Fig. For each spin in the first and the last 
columns, the effect of these dynamic boundary conditions is to effectively decouple this spin from its up neighbor. In 
this way we mimic a concave, stair-like border (as the one found by the domain wall in the circular nanoparticle at 
some stages of its evolution, see Fig. [T5|) . with a fixed distance between both concave borders. 

We initialize the system with a flat domain wall between the metastable and stable phases. Under the action of the 
magnetic field, the domain wall propagates toward the metastable phase. Due to the concave boundary conditions in 
the X direction, the interface propagates faster near the boundary. After a short transient, the initially flat domain 
wall reaches an stationary state, with an almost constant (up to small fluctuations) non-zero curvature, which depends 
linearly on the size Lx, see Fig. In order to emulate an infinite system in the ij direction, we shift the observation 
window whenever the domain wall gets close to the system extreme in the propagation direction, always keeping the 
interface inside the system. In practice we generate a new region with up (metastable) spins in one extreme of the 
system whenever the above condition is met, eliminating an equivalent region of down (stable) spins in the opposite 
extreme. (In fact, in order to give time for the newly introduced spin region to relax to the typical state of the 
metastable phase at the given T, p and we perform the shift well in advance before the domain wall reaches the 
boundary.). 

The above defined system evolves via avalanches, whose distribution we can measure. In this way we obtain a 
reliable measure of avalanche statistics for a domain wall with (approximately) constant curvature, see the right 
graph in Fig. 1161 In particular, we observe that domain walls with constant curvature exhibit avalanches with a 
well-defined typical scale. 
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